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Abstract. We present long-term observations of 
several VHE (E > 100 GeV) blazars with VERITAS, 
together with contemporaneous Swift and RXTE X- 
ray data. The observed targets include Mrk421, 
Mrk501, 1ES 2344+514. Strong flux and spectral 
variability is seen in Mrk421 on nightly time-scales 
between January and June 2008, revealing a highly 
correlated X-ray to GeV/TeV connection. Modest X- 
ray variability is evident in Mrk501. Observations 
of 1ES 2344+514 in December 2007 show VHE 7-ray 
and X-ray flux doubling on nightly time-scales. 

Keywords: Gamma-ray Astronomy, Active Galactic 
Nuclei, VERITAS 

I. Introduction 

Blazars (BL Lac objects and Flat Spectrum Radio 
Quasars) are active galactic nuclei (AGN) with a rela- 
tivistic plasma jet oriented close to the line of sight Q. 
These objects exhibit rapid variability and have broad- 
band spectral energy distributions (SEDs) characterized, 
in a v¥ v representation, by a synchrotron component 
extending from radio to X-ray frequencies, and a second 
component peaking at 7-ray frequencies due to either 
inverse-Compton radiation or from hadronic processes. 
BL Lac type blazars are further sub-divided based on 
the peak frequency of the synchrotron emission as 
low-, intermediate-, or high-frequency-peaked BL Lacs 
(LBLs, IBLs, and HBLs). Currently, 20 HBLs from a 
total of ^25 VHE blazars are detected at very high 
energy (VHE, E >100 GeV) 7-raysQ Elsewhere in these 
proceedings, the VERITAS blazar observing program 
is described Q, and the discovery of VHE emis- 
sion with VERITAS from two HBLs (RGB J0710+591 
and 1ES 0806+524), and from two IBLs (3C66A and 
W Comae) is highlighted |3|. 

Multiwavelength campaigns on BL Lacs are important 
to sample the broadband flux and spectral variability 
on time-scales ranging from minutes to months. Recent 
joint observing campaigns with VERITAS on two IBLs 
are presented elsewhere in these proceedings. Light 
curves and the broadband SED of 3C 66A are presented 
from contemporaneous VERITAS, Fermi, X-ray, and 
optical data between September and November 2008 [4J. 
Observations of W Comae in June 2008 reveal X-ray 
and VHE 7-ray flaring, with a marginal detection above 
100 MeV from AGILE observations J5). In this paper 

'TeVCat catalog of VHE 7-ray sources: http://tevcat.uchicago.edu 
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Fig. I. Time-averaged VERITAS photon spectra of Mrk421 for 
discrete flux levels. A power law with exponential cutoff model 
dN/dE = I -(E/lTeV) r -exp(-E/E cut ) is fit to each spectrum. 



detailed multiwavelength results on the bright and well- 
studied HBLs: Mrk421, Mrk501, and 1ES 2344+514 
are presented. Multiwavelength results on the recently 
observed HBL RGB J07 10+591 will also be presented 
at the conference. 

II. VERITAS and X-ray Data Analysis 

VERITAS is an array of four imaging atmospheric - 
Cherenkov telescopes located in southern Arizona. The 
array is sensitive over the energy range of 100 GeV 
to >30 TeV, and can detect (5 a level) a source flux 
of 5% of the Crab Nebula flux in ~2.5 hours. All 
observations presented here pass run quality selection 
criteria, which removes data taken during poor weather 
conditions or non-standard hardware configurations. The 
wobble mode, where the source is positioned at a fixed 
offset of 0.5° from the camera center, was used for 
all observations presented here. Standard data reduction 
and 7-ray selection cuts are applied [6]. The results 
presented here agree well with those performed using 
an independent VERITAS analysis package. 

Long-term X-ray observations with RXTE PCA Q 
and Swift XRT [8] were taken contemporaneously with 
the VERITAS data. All RXTE PCA data were taken 
with only PCU2 operational. For all Swift XRT data on 
Mrk421 the observations were taken in window timing 
(WT) mode, while for 1ES 2344+514 the XRT data were 
all taken in photon counting (PC) mode. Data reduction 
is performed with the HEAsoft 6.5 package, following 
the standard methods 0, iflOl . In particular, annular 
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Fig. 2. VHE 7-ray and X-ray light curve of Mrk421. Shown in the top panel are VERITAS F(>300 GeV) fluxes from observations of 20 
minutes each. The middle panel shows the 2-10 keV fluxes from observations with RXTE PCA (circles) and Swift XRT (open squares). The 
bottom panel shows the X-ray spectral indices (a parameter, see text) from log-parabolic fits to the 0.4-10 keV Swift XRT spectra and 3-20 keV 
RXTE PCA spectra. 
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Fig. 3. VERITAS flux F(>300 GeV) from 20 min exposures versus 
the X-ray flux for observations of Mrk421 taken within ±5 hours. 

source regions are used to extract the Swift XRT data 
in PC mode when photon pileup is evident [9|. 

IE. MRK421 

Mrk421 is a nearby (z = 0.031) HBL and is highly 
variable at UV to VHE 7-ray energies. VERITAS ob- 
served Mrk421 between January and June 2008 for 
a total live-time of 43.6 hours. An average analysis 
energy threshold of 260 GeV is determined for the 
observations that span a range of zenith angles from 



6°-56°. In the whole data set an excess of ^30000 7- 
ray events is detected, corresponding to a significance 
of >270 a. Previously observed VHE 7-ray spectral 
variability correlated with flux level ifTTI is investigated 
by dividing the VERITAS data into subsets according 
to flux level. Figure Q] shows the VERITAS spectra 
in six subsets with best-fit curves for a power law 
with exponential cutoff model of the form dN/dE = 
I • (E/lTeV)~ r • exp(-E/E cut ). The nightly X-ray 
spectra from RXTE PCA and Swift XRT data are best 
fit with an absorbed log-parabolic model, which uses 
a fixed column density iflOl and an energy dependent 
photon index T = a + b- Log(E/E Q ). The mean reduced 
X 2 values from log-parabolic fits to the RXTE PCA 
and Swift XRT spectra are 0.87 and 1.20 compared to 
reduced \ 2 values of 1.35 and 1.80 for an absorbed 
power law model. 

Figure [2] shows on the top panel the nightly VHE 
7-ray fluxes F(>300 GeV) from nightly VERITAS ob- 
servations. Shown on the middle and bottom panels are 
X-ray fluxes F(2-10 keV) and spectral indices from 
RXTE PCA and Swift XRT observations. Strong flux 
variability on nightly time-scales is seen over the entire 
six month period in 2008, with exceptionally bright X- 
ray flaring in March to April, and VHE 7-ray flaring 
in early May. A measure of the integrated level of 
flux variability is the fractional root-mean-square (rms) 
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Fig. 4. VHE 7-ray and X-ray light curve of Mrk501. Shown in the top panel are the nightly VERITAS F(>300 GeV) fluxes. The middle 
and bottom panels shows the RXTE PCA 2-10 keV fluxes and power law indices, respectively. 



variability amplitude F var IIT21 . Significant VHE 7-ray 
flux variability of F var = (66.5 ± 0.6)% and X-ray 
flux variability F var = (57.4 ± 0.3)% was observed. 
Figure [3] shows the VHE 7-ray flux versus X-ray flux 
for observations taken within ±5 hours of each other. 
The highest near-simultaneous X-ray fluxes are from the 
March to April period. Unfortunately, when VERITAS 
measured the highest VHE 7-ray fluxes on May 2 and 
May 3 there were no simultaneous RXTE PCA or Swift 
data. Combining all VHE 7-ray flux versus X-ray flux 
points, the correlation coefficient is r = 0.62 ± 0.07. 
Clear X-ray spectral hardening with increasing 2-10 keV 
flux is shown in figure 12 with the log-parabolic index 
parameter a ranging from ^ 1.8-2. 8. 

IV. MRK501 

The nearby HBL Mrk501 (z = 0.034) was discovered 
at VHE 7-ray energies by the Whipple 10 m telescope 
in 1996 HI. VERITAS observed Mrk501 between 
April and June 2008 for a total live-time of 6.2 hours. 
An excess with the statistical significance of 22 a is 
measured for the total data set. Figure [4] shows the 
nightly VERITAS VHE fluxes F(>300 GeV) and RXTE 
PCA 2-10 keV fluxes and power law indices. Marginal 
variability is seen in the VERITAS F(>300 GeV) fluxes 
with F var = (14 ± 12)%, while moderate 2-10 keV 
flux variability is evident with F var = (12.7 ± 0.4)%. 
More detailed results from this campaign are presented 
elsewhere in these proceedings lfl4l . 



V. 1ES 2344+514 

1ES 2344+5 14 is another close HBL (z = 0.044), and 
was first detected at VHE 7-ray energies by the Whipple 
10 m telescope in 1995 1151 . Figure [5] shows the nightly 
VHE 7-ray and X-ray light curve of 1ES 2344+514 from 
VERITAS, RXTE PCA, and Swift XRT observations. A 
strong VHE 7-ray flare is seen on December 7, 2007 
(54441.12 MJD) at a flux F(>300 GeV) correspond- 
ing to 48% of the Crab Nebula flux. The measured 
increase in flux of a factor of 1.92 ± 0.53 between the 
previous night and the flare night shows the first clear 
evidence of ~day time-scale VHE 7-ray variability from 
1ES 2344+514 since the initial Whipple 10 m detection 
in 1995. Excluding the December 7 flaring event, the 
average F(E>300 GeV) is 7.6% of the Crab Nebula flux. 
For the full VERITAS data set a high level of variability 
F V ar = (75 ± 10)% is implied. Excluding the flare night, 
a F var = (34 ± 16)% is determined. 

Figure [5] (lower panels) shows the 2-10 keV flux and 
photon index T measured over 3-20 keV from RXTE 
PCA and 0.4-10 keV from Swift XRT data. The X- 
ray flux is shown to be highly variable throughout the 
campaign, with F var = (51 ± 1)%. In December 2007, 
large amplitude flaring is evident with flux doubling 
time-scales of -1 day. A 2-10 keV flux of (6.58 ± 
0.37) x 10~ n erg cm' 2 s _1 is seen from the Swift 
XRT data on December 8, 2007, representing the highest 
X-ray flux ever measured for 1ES 2344+514. Figure [6] 
shows the VERITAS flux F(>300 GeV) versus RXTE 
PCA and Swift XRT 2-10 keV fluxes for nights with 
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Fig. 5. VHE 7-ray and X-ray light curve of 1ES 2344+514. Shown in the top panel are the nightly VERITAS F(>300 GeV) fluxes. The 
middle panel shows the 2-10 keV fluxes from observations with RXTE PCA (circles) and Swift XRT (open squares). The bottom panel shows 
the power law indices T from the 0.4-10 keV Swift XRT spectra and 3-20 keV RXTE PCA spectra. 
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Fig. 6. VERITAS 7-ray flux F(>300 GeV) versus X-ray 2-10 keV 
flux for 1ES 2344+514 observations from nights with RXTE PCA 
(circles) and Swift XRT (open squares) data. 



observations in both energy bands. A Pearson coefficient 
of r = 0.60 ±0.11 is calculated for the VHE 7-ray to 
X-ray flux points, suggestive of correlated variability. 
Further results and discussion on this campaign are in 
press 1 16]. 

VI. Conclusions 

Joint VHE 7-ray and X-ray observing campaigns on 
the bright HBLs Mrk421, Mrk501, and 1ES 2344+514 
with VERITAS reveal significant flux variability on 



nightly time-scales. X-ray spectral hardening is shown 
at increasing flux levels. Results from synchrotron self- 
Compton (SSC) modeling of the broadband SED will be 
presented at the conference. 
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